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Abstract: We investigate the spacetime-dependent condensation of the tachyon in effec- 
tive field theories. Previous work identified singularities in the field which appear in finite 
time: infinite gradients at the kinks, and (in the eikonal approximation) caustics near lo- 
cal minima. By performing a perturbation analysis, and with numerical simulations, we 
demonstrate and explain key features of the condensation process: perturbations generi- 
cally freeze, and minima develop singular second derivatives in finite time (caustics). This 
last has previously been understood in terms of the eikonal approximation to the dynam- 
ics. We show explicitly from the field equations how this approximation emerges, and how 
the caustics develop, both in the DBI and BSFT effective actions. We also investigate the 
equation of state parameter of tachyon matter showing that it is small, but generically non- 
zero. The energy density tends to infinity near field minima with a charateristic profile. 
A proposal to regulate infinities by modifying the effective action is also studied. We find 
that although the infinities at the kinks are successfully regularised in the time-dependent 
case, caustics still present. 
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1. Introduction 



In Type II string theories, unstable branes of dimension p decaying into the stable ones 
with lower dimensions is well described by the Dirac-Born-Infeld (DBI) type effective action 

S = - J dP+'xV{T)./TT^, (1.1) 

where y = i']^'^ d^Td^T , V{T) is the asymptotic potential, and we use units in which 
a' = 1. At the beginning of the condensation, the system has zero tachyon field T = 
and is located on the top of the potential V{T). It is unstable and will roll down when 
driven by a small perturbation. The field evolves towards the minimum of the potential 
y = 0, and the condensation ends. The final state corresponds to the disappearance of the 
unstable branes, replaced by "tachyon matter" whose properties are not well-understood. 

The condensation process in the inhomogeneous case has also been studied. It has 
been realised that the equation of motion from the DBI action leads to solitonic solution 
with kinks and anti-kinks [@, ^, ^ 01- At the kinks and anti-kinks, the field remains zero, 
and become daughter branes of one dimension lower at the end of the condensation. While 
the field in between them grows with time and tends to infinity. The huge difference of the 
field values between these two areas makes the field gradient at kinks and anti-kinks very 
large and in fact can be semi-analytically shown to lead to infinity in finite time 0]. The 
existence of this kind of singularity prevents straightforward numerical integration of the 
space-time dependent field equation IC]. 



In between the kinks and anti-kinks where the field is approximately homogeneous, an 
eikonal approximation 1 + y = reveals that the second and higher derivatives of the field 
can reach infinity in finite time. This is interpreted as the formation of caustics |^. 

The potential application of the tachyon condensation in the inflation scenario has 
been investigated, regardless of the various problems |11, 12, 13, 14, 15, 16]. In the tachyon 
inflation model, agreement with the observational data can be achieved [17|. A complex 
tachyon field can produce cosmic strings after infiation |18, 19, The appearance of 
the pressureless tachyon matter at late stage of the condensation [^, ^ is also a feature 
of this model. The fiuctuations of the tachyon matter in both the free and the expanding 
cosmology background [22, |l^, |l^, 17] have been studied, and shown to increase linearly 



with time, just like ordinary pressureless matter. 

In this paper, we study the condensationg process using the field equations, demon- 
strating and explaining, both analytically and numerically, key features of the inhomo- 
geneous tachyon condensation process. Perturbations around the homogeneous solutions 
T = t "freeze", that is, show little change after a short initial relaxation period. What 
change there is reduces the curvature at maxima, and increases it at minima, until it be- 
comes singular in finite time. This is interpreted as the formation of a caustic in the eikonal 
solution. 

We also show explicitly from the field equations, backed up with numerical simulations, 
that the quantity 1 + y relaxes exponentially to zero, thus accounting for the accuracy of 
the eikonal equation. 
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We also investigate the energy density and pressure of condensing tachyon matter, 
showing that it is not quite pressureless, and that the energy density diverges in a charac- 
teristic l/|Aa;| from near the developing caustic. 

One can ask if the instabihty is an artifact of the DBI effective action. To give a partial 
answer, we also study the boundary string field theory (BSFT) effective action 
showing that this approaches the eikonal equation eventuahy. 

Finally, we show that a modification of the action proposed by ||5| to regulate the 
gradient at a static kink does prevent the gradient of the field reaching infinity at the kinks 
in the time-dependent case, but does not prevent caustic formation between kinks. 

The paper is constructed as follows. In Sec. 2, we give the equation of motion and its 
solutions in simple cases. In Sec. 3, we present the 1 + 1 dimensional solutions near kinks 
and near extrema respectively. In Sec. 4 and Sec. 5, we discuss the dynamics approaching 
the vacuum respectively in the DBI and the BSFT effective theories. Based on the results, 
the features of the tachyon matter towards the end of the condensation are investigated in 
Sec. 6. In Sec. 7, we discuss the generalised DBI action. 



2. Equation of motion 

From the DBI action (|1.1| ), we can write down the equation of motion: 

□T-y) (l + y) = Vra^(l + y), 
where □ = r]^jjd^d^ . An equivalent expression to the equation of motion is 
T = f I 2TViTV'T + (1 + y)V^T - ViTVjTV'V^T -y{l + y) 



(2.1) 



(2.2) 



where / = 1/(1 + ViTV^T), = did^ and i, j = 1, ■ ■ ■ ,p. V{T) is the runaway tachyon 
potential and its field derivative is V = dV{T)/dT. A suitable choice of potential, derived 
from the boundary conformal field theory on the worldsheet [ 25 , |2^ , is 



V 



cosh (pT) ' 



(2.3) 



where the constant /3 = 1 for the bosonic string and /? = l/y/2 for superstring. 
The energy-momentum tensor is: 



V{T) 



Vl + dT- dT 



Vl + dT- dT 



(2.4) 



2.1 The static case 



For the 1-dimensional static case, the solution to the field equation is most easily found by 
noting that conservation of energy- momentum requires that the pressure is constant: 



-V 



'11 



Vl + T'^ 



-Vn. 



(2.5) 
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where Vq = V(Tq) is the minimum potential when the maximum field Tq is achieved at 
X = xq. The maximum field satisfies the condition T'\t=To = 0- 

For the inverse hyperbolic potential ( |2.3D , the equation is solvable: 



T(x) = — sinh ^ 



— — 1 sm{PAx) 



(2.6) 



where r = Vo/Vm and Ax = x — Xm- It represents an array of solitonic kinks and anti-kinks. 
Its gradient is: 



T' = ±4 



r2 + tan^ (/3Ax) ' 



(2.7) 



For fixed Vm, T' tends to infinity at kinks and anti-kinks, where Arc = nvr//?, as Vq — > 0, 
which is consistent with the slow motion analysis on the moduli space of W\. 



2.2 The homogeneous case 

In this case, the energy density is conserved and constant: 

V 



too 



Vl-T2 



which restricts |T| < 1. 

For the potential (|2.3|), the time dependent solution is: 



and its time derivative is: 



isinh-i[^^ - lcosh(/3t)], {E < Vm), 
isinh-i[Jl-^sinh(/5t)], {E>Vm), 



± 
± 



1-/2 



coth2(/3t)-Z2 • 



{E > Vm), 



(2.8) 



(2. 



where / = E/Vm- For E > Vm, T = 0, \T\ = y^l - l/P at t = and |r| 
t — > oo. 



(2.10) 



oo, \T\ — > 1 as 



3. Solutions near kinks and extrema 

In this section, we study the behaviour of solutions near stationary kinks and extrema, in 
order to gain insight into the singular dynamics, and to make contact with previous work 
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3.1 Around kinks 

Following P], we expand the field around a stationary kink or anti-kink located at x = Xm'- 



T{t,x) = a{t)Ax + -c{t)Ax^ + —e{t)Ax^ + 
6 120 



(3.1) 



where Ax = x — Xm- 

For the potential (^), V /V = -p'^T. For the potential V /V ~ -/S'^T when 

|r| is small. We use this relation and the expansion (3.1) in the field equation around kinks 
and anti-kinks, where |T| is small. The comparison of coefficients gives the equation: 



2 2aa + c 
a = 13 a + — — — 
1 -I- 



(3.2) 



We are more interested in the solution at late time when a{t) becomes large. If we can 
neglect the contribution from c(t), 



aa = (3^a^ + 2d^. 
To solve the equation, we set a = 1/z, giving 



with solution 



z = z\ COS (/3t) -|- Z2 sin (/?t). 
By a suitable choice of the time coordinate, the solution of a can be written 

ao 



a(t) 



cos (/3t) 



(3.3) 

(3.4) 
(3.5) 

(3.6) 



where ao is constant. If we set = 7r/(2/3), we approximately have: a ~ l/(t* — t\ which 
is consistent with the result of a(t) grows to infinity in a finite time. 

More precisely, we now also consider the evolution of the coefficient c(t). With the 
same approximation of a being large, we can obtain the equation: 



a^c = 8adc -|- (7/3^a^ — 6ad)c — 6/3^ ad^ -|- e. 
Combining with Eq. ( |3.3| ) and neglecting e(t), we find 

c(t) = -13 a = -— . 

cos [pt) 

Therefore, the solution near the kinks and anti-kinks can be expressed as: 



T{t,x) = a{t) 



Ax - -d^Ax^ + 



ao sin {13 Ax) 
^3 cos {(3t) 



(3.7) 



(3.^ 



(3.9) 



This has the following properties at x = Xm- T = 0, T = 0, T'{t oo) — > oo, T = and 
T" = 0, which are consistent with the discussion in Sec. 2. 
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3.2 Around extrema 

The field around a stationary peak or trough at x = xq satisfies the condition T'\x=xo = 0. 
The expansion around xq is written: 

T{t, x) = To{t) + h{t)Ax^ + ^d{t)Ax^ + ■■■ , (3.10) 

with Ax = x — xq. For convenience, we only consider the behaviour around the field 
extrema whose values are positive To(t) > 0. 

For the potential ( ^.31) , V' /V ~ — /? when T is positive. Inserting the expansion form 
in the field equation (12.21) , we get: 

Tb = (/3 + 6)(l-T2). (3.11) 

The final value of Tq is decided by the relation between b{t) and /?. Here we only consider 
two special cases. When b = — /3, Tq is constant and it corresponds to the static case. This 
is consistent to the static analytical solution ( |2.6[ ), which has T" = — /3 at peaks. Secondly, 
if |6| <^ P, there is an explicit solution for Tq 

To = tanh(/3t). (3.12) 

The solution satisfies Tq £ [0,1]. Otherwise, the solution can be formally expressed as: 
Tq ~ 1 — {A/2)e~'^^^~^^^^ , where j4 is a small positive value. 
The equation for the coefficient b{t) is: 

b = 2(3b^ - 2(/3 - b)tQb. (3.13) 

Since we are considering the peaks on the positive field side Tq > 0, Tq should be positive 
and tends to 1. Making use of the approximation Tq ~ 1 for late time, Eq. ( |3.13| ) can be 
rewritten as: 

dt{b-b^) = -2l3{b-b^). (3.14) 

Hence, 

6 = 62 + (^gxp(-2/3t), (3.15) 

where C = bQ — b"^ is a. constant, with 69 = b{t = 0) and Jjq = {db/ dt)\t=Q. We have made 
numerical solutions of ( |3.15| ). The distribution of the final values as a function of 69 and 
60 is presented in Fig. ||. From these results, we see that a peak (60 < 0) is most likely to 
flatten with time, but a trough (60 > 1) can deepen and become singular, and will always 
do so when 60 is positive. There is a sign of the developement of a caustic. 

Focus on the late time behaviour of b{t). In terms of Eq. ( 3.14| ), we will have ast 00 



b = b\ (3.16) 

So we can have the solution of b 

m = (3.17) 

where 60 = b{t = 0). Thus b blows up at t = Iq when 60 > while 6 — > when 60 < 0. But 
b > for both cases. This is consistent with the numerical solutions prsented in Fig. (||). 
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b(t) 



Figure 1: The time evolving curves of b{t) vs b{t) with different initial values {bo, bo) in regions 
where b > ~f3. The blue, solid lines indicate that b{t) at large t and the red, dashed lines 
indicate that b{t) tends to infinity. 

3.3 Caustics around extrema 

In terms of the relations in Eq. 
late time as: 

l + y^l-t^-{b- b^)£\x^ ~ - C £\x^)e-'^<^\ (3.18) 

If 6 > —{3, we see that 1 + y ^ 0, first noted in |^. If we are in the region 6 < 6^, C is 
negative and 1 + y approaches zero from above. 

Based on this relation, the authors of Q argue that there should be caustics forming 
around the peaks or troughs. They rewrite the relation ( p.l8| ) as an eikonal equation: 

52 - 5'2 = 1, (3.19) 



3.12) and (3.15|), we can estimate the quantity 1 + y at 



with the approximation S* ~ T and develop a series of characteristic equations for ( 3.19| ) 
to give 

(l+S"2(0,x))3/2^ 

Thus S" blows up in a finite time i = (1 + 5'2(0, x)f/'^/S"{{), x) when 5"(0, x) > 0, which 
is interpreted as the appearance of caustics or regions where the tachyon field becomes 
multivalued. We also note that when S"(Q,x) < there are no caustics, with S"{t,x) 
as t — > 00. 

Our results are consistent with the HJ analysis above: caustics form near troughs 
{b > 0) in a time 6q ^, where 60 — S"{0,x) at a trough, and peaks tend to flatten (b 0). 
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4. Perturbations 



In this section, we will consider the evolution of perturbations around a condensing homo- 
geneous tachyon field in an arbitrary dimensional Dp-brane. This will also approximately 
describes the evolution of the field between well-separated kinks and anti-kinks. A general 
form is: 

Tit,x') = Toit)+T{t,x'), (4.1) 

where T{t,x'^) is the small perturbation along all spatial directions along the Dp-brane 
with z = 1, • • • ,p. 

A trivial case is: To{t) = 0. Around To{t) = 0, -V /V ~ (3^t for both potentials 
( |2.3| ) and (|5.lD . Inserting T{t,x'^) = T{t,x^) in the field equation ( |2.2[) , we can get the 
perturbation equation: 

VV - f + /3V = 0, (4.2) 
as expected for a tachyon field with mass squared — 

4.1 Linear perturbations around a condensing field 



Plugging the ansatz (4J.) in the field equation with the hyperbolic potential (|2.3| ), we have 
the equation of r 

[To - /3(1 - t^)] + [f + 2(3tof - (1 - To') VV]+, (4.3) 
[Pf^ - (/? - ro)Vr • Vr - 2To(Vt • Vf - fVV)] 

Solving the zeroth order equation covered in the first square bracket, we get the solution 
To = tanh(/3t), as given in Eq. ( 3.121) . 

We first consider the linear perturbation covered in the second square bracket. This 
situation applies if all the first order terms remains dominated compared to the second 
order terms, though this is always not true as we can see later. 

The solution of r in the second square bracket is separate, and writing r(t, x) = 
f{t)g{x), and we find: 



V^g + k'^g = 
/ + 2/3tanh(/3i)/ + k^sech^ifSt)/ = 



where /c^ = /cj/c* and fc* are arbitrary constants. 
By making the reparameterisation 



^ (^<P<1), (4.5) 



we can rewrite the second equation in Eq. (4.4) as 



P(1-P)0 + ^V = O. (4.6) 
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where and Ki = ki/(3. It is a typical hyper geometric differential equation. The 

equation has two linearly independent solutions, the simpler one of which is given as |27| 



f{t)=p{t)2Fi ^ , ^ ■2;p{t) , (4.7) 



where is the hypergeometric function. 

We will focus on the late time behaviour of r when Tq is closed to 1. In this situation, 
the second equation in Eq. (p^ ) becomes / + 2/3/ + 4/c2 exp (— 2/3t)/ = 0. With the 
reparameterization = exp (— 2/3t), we have: 



v-j:{ + ^^f = o. (4.8) 



The equation is a Bessel style differential equation . Its solution is given as a conbination 
of two Bessel functions. The whole solution at late time can be expressed as: 

r flPh r 1 



. (2vr)P - 

where ri(/cj), r2(/cj) are constant parameters corresponding to the mode ki and 

/i(^'2)(t,/fc) = e-f''H[''^\2Ke-'''). (4.10) 

' [s] are the Hankel function, which are the linear combinations of the Bessel function 
of the first kind Ji(s) and the second kind Yi{s): h[^\s) = Ji(s) + ^^1(5) and h[^\s) = 
Ji(s) - iYi{s). They satisfy h[^\s) = h[^\s)* . In order to get r = t* , Ti{ki) = T2{ki)* . 
Then we can set Ti{ki) = To{ki)e'^^^^^ and T2{ki) = To{ki)e~'^^'^^\ 

Here, we consider the real tachyon field T. The perturbation ( [4.S| ) should be real as 
well and can be expressed as, denoting s = 2Ke~^*: 

r fpu o 

T{t,x') = J J^^Mh) [Jiis)cosihx' + 6iki))-Yiis)sm{hx' + 6iki))] . (4.11) 

For small parameters < s ^ \/2, 

Ji(s)~^, Y,{s)c^- — . (4.12) 

4 TVS 



Thus, at late time, Eq. ( |4.11 ) approximately reduces to: 



The result indicates that the perturbation corresponding to each mode ki freezes at late 
time and that shorter wavelength modes are damped more than longer wavelength. 
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4.2 Comparison with the numerical simulations 

We make numerical simulations of the perturbations in 1 + 1 dimensions based on the 



Hamiltonian formalism [29, 3C, y]. Define the momentum conjugate to T 



U='4 = ^. (4.14) 



Then the equations of motion can be expressed as 



T=^^v^^. (4.16) 

Vn2 + y2 ^ ' 

The simulation is implemented on a spatial lattice of 1280 points. The lattice spacing 
and the timestep are set to be respectively: 5x = 0.05 and 5t = 0.01. We use the symmetric 
difference for the first order field derivatives and a 3-point stencil for the second order field 
derivatives. We adopt the second order Runge-Kutta method for the time update. In the 
simulation, the constant (5 = l/\/2- 

To get the perturbation, we first need to produce two sets of simulation data with two 
different initial conditions respectively: T[t = 0, x) = D and T{t = 0, x) = D + 5{x), where 
D is the field value at t = 0: D = ro(t = 0), and 5{x) is the random initial perturbation: 
5{x) = T{t = Q,x). The difference between them gives the 1 + 1 dimensional evolution 
surface of r(t, x). In Fig. ^, we present the plots of the perturbation at t = 8,24,38.4 
respectively. We can see that the perturbation almost "freezes", especially at peaks and 
troughs, which is consistent with the linear perturbation analysis in the previous subsection. 
However, the second derivative V'^t can be seen to decrease in magnitude near peaks, and 
increase near troughs, until the simulation breaks down (signalled by a spike in r at around 
X ~ 27 and x ~ 53). These features are consistent with the analysis in Section 3, but not 
with the predictions from the linear perturbation equation. Hence, we need to take into 
account higher order perturbation terms to understand this behaviour. 

4.3 Non-linear perturbations and development of caustics 

We start with a discussion of the behaviour of 1 + y. For large T, we can assume that 
T — > 1, T ^ and V /V ~ —j3. The equation of motion ( p.l| ) becomes 

-(1 + y) ~ -2(v2r + (3){i + y) + 2V,rVjrv,v,r - t—{v^t). (4.i7) 

If the last two terms are negligible compared to the first one, then 1 + y should decrease 
exponentially with time to zero. This accounts for the numerical observation first made in 

To check the expectation that 1 + y ~ 0, we present the plots of 1 + y in Fig. (|3|) for 
a field T(t, x) = D + t + T{t,x), where the plots of the perturbation r are given in Fig. 
The result indicates that 1 + y really decreases with time to very small positive values 
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>< 




Figure 2: The numerical plots of the perturbation T{t,x) based on Eqs. ( 4.15| ) and (4. If). The 
initial homogeneous field is set to be D = 1. 




Figure 3: The plots of log(l + y) for the field T{t, x) — D + t + T{t, x), where D = 1 and the plots 
of T{t, x) are given in Fig. (||). 



exponentially in the early stage. However, the simulation can not continue at late time 
due to instabilities emerging at the positions where r eventually gets discontinous in Fig. 
@, i.e., near troughs. This implies that the instability should correlate with the caustics 
formation. In what follows, we will determine how it happens in the perturbation method 
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Including the quadratic order terms of r in the second square bracket in Eq. ( [4.31 ), we 
have in the hmit Tq ^ 1 

dt[2f - [Vrf] = -(3[2f - (Vt)^] +f- f(2VV + (3f). (4.18) 



It is consistent with Eq. ( 3.14| ) if setting r = b{t)Ax'^/2 in p = 1 case. When the last two 



terms on the right hand side of the equation are less important than the first one, we will 
get at late time 

l + y~2f-(VT)2 = 0. (4.19) 

Specially, f = when VjT = 0, i.e., at exactly the peaks or troughs, the perturbation 
"freezes". But away the peaks and troughs, the perturbation increases with a increasing 
rate equal to (Vr)^/2. 

Let us see what features we can have from the equation f = (Vt)^/2 around peaks 
and troughs respectively. Since r "freezes" at exactly the peaks and troughs and grows 
away from them, the field gradient VjT decreases near peaks and increases near troughs. 
Therefore, the curves of r tend to flatten near peaks untill V^r = 0, where f = and so 
r stops growing. While near troughs, the curves of r sharpen, with |Vjr| increasing. This 
accelerates the growth of r. In this case, r, V^r and V^r all will tend to infinity, which 
accounts for the caustic formation and so the instability near the vacuum. 

5. The BSFT action 



The effective action derived from BSFT is given by [23, 24]: £ = —V{T)F{y) with the 
potential 



F(T) = exp(^-^J, (5.1) 
and the kinetic part 

Fly) = 1^01 = 2^.-VB(,), (5.2) 

where y = d^Td^T as before and B{y) is the beta function. It is easy to see that F >{) 
when y > —1/2. Aspects of the dynamics arising from the the BSFT action have been 
given in |3^ . 

Using r(l + u) = nr(n), we can derive the recursion relation 

B{y) = 2^- /^^^"^,^^^,^;"^"^!! B{y + m), m = 1, 2, 3, • • • (5.3) 
(y + m - l)(y + m - 2) • • • (y + l)y 



From Stirling's approximation, B{u) ~ 2^ ^" y^vr/ti for large u. Then the BSFT action can 
be expressed as 

(y + ^)(y + ^)---(2/ + ^) 



F{y)^i^ ' ; ^ , x W— — . (5-4) 

[y + m - l){y + m - 2) ■ ■ ■ [y + 1) \ y + m 



and we will find it useful to take m ^ —y. 
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For m = 1, we approximately have 

F{y) ~ + (5.5) 

when 1 + y ^ 1. This is the DBI action. So the behaviour of the BSFT action ( ^.21 ) for 
large 1 + y should be similar to that discussed in previous sections. In what follows, we 
will discuss the behaviour near the vacuum when y is small or negative. 

5.1 Time evolution 

We first explore the homogeneous approach to the vacuum. The equation of motion in this 
case is ^ 

F-2yF' = -, (5.6) 

where y = — T^, F' = dF/dy and E is a constant. 

The condensation starts with a small velocity |T| . Let us consider the case when T = 0. 
With the relation of the gamma function T{2u) = 22«"ir(u)r(u + l/2)/V¥, we can rewrite 
the kinetic part F{y) as 

F(rt = (5.7) 

r(y+ 2) 

Thus, F{0) = 1 since r(l/2) = /tF. This gives V{T) =Eatf = Oin terms of the equation 
of motion. 

Now we need to find where the tachyon potential gets the minimum value V(T) = 0, 
which corresponds to the vacuum. Since y = — < in the homogeneous case, we can 
change the action to the form by using r(l — u)T(u) = vr/ sin(7rn) 

F{y) = -22s/+i^cot(7ry)— ^ ~ _^cot(^y)^. (5.8) 
Then the homogeneous equation of motion ( |5.6| ) becomes 

(5.9) 



2(-^j/)^ ^ E 
sin^(7ry) V 



It is easily seen that V = at y = —n, where n = 1, 2, 3, • • • . The behaviour in taking the 
limit — > n is, 

nr)~fyj(r2-n)^ (5.10) 

The equation implies that the tachyon potential vanishes for either ^ n or n — > 00. We 
will assume that only the first pole (y = — 1) in Eq. ( ^.91) is relevant, as the condensation 
starts with y ~ 0. 

5.2 Dynamics near the vacuum 

The full equation of motion for a general F{y) is 

(^2^d,Td^T + r^,^y^d'^T+^l^y-\^^=0. (5.11) 
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Note that V = dV/dT, F' = dF/dy and F" = d'^Fjdy^. 

We now consider the dynamics approaching the vacuum when y — > — 1 from above. 
Since FjF' = l/(lnF)' and F"/F' = (lnF)'7(lnF)' + (lnF)', we only need to know (InF)' 
and (InF)" for deriving the equation of motion. In the hmit y — > —1, we approximately 
have for any m >2 from Eq. (| 



(InF)'---^, (InF)"-— (5.12) 

The corresponding expressions for the DBI action are (InF)' = 1/[2(1 + y)] and (InF)" = 
— 1/[2(1 + y)^]. Inserting them in the equation of motion, we have 

□T-^^ (l + y)~25^Ta^(l + 2/), (5.13) 

So it is very similar to Eq. ( p.l| ) for the DBI action. The only difference is a change of the 
coefficient on the right hand side. Therefore, there should be similar consequences for the 
dynamics approaching the vacuum from the BSFT action to those from the DBI action. 
That is, the perturbations "freeze" and the second derivative becomes discontinuous. 



6. Tachyon matter 

Reverting to the DBI action, the energy density and the pressure are given in flat spacetime 
by Eq. (0): 

p = roo = ^f=^y(r), (6.1) 

l + y^-^V^ (ao^ 
P. = Tu = -V{T) ^^^^ = g + Y^^P, (6.2) 

where there is no sum on i and q = C = —V{T)yJ\ + y, and so 

_Pi _ 1 + ^ , (V»^)' 

p~ i + (vr)2 + i + (vT)2- ^^-"^^ 

From the above formulas, it is easy to see that the system at the beginning of the 
tachyon condensation has positive energy density and negative pressure everywhere: pi = 
—p = —Vm, the characteristic of the interior of a brane. After the condensation starts, the 
behavior in different areas becomes different. 

Consider the case that the decay of the Dp-brane happens only in the direction to 
produce a D{p — l)-brane. Then the decay process can be described by a kink-anti-kink 
tachyon solution along the direction. 

First, we consider the appearance of a D(p — l)-brane, which we locally choose to be 
orthogonal to the x^-direction. Since T = and |ViT| oo towards the end of the decay 
at kinks, the energy density near the kinks px — > oo. If the approximate solution found 
in Sec. 3.1 is to be believed, this divergence happens in finite time, in accordance with 
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the analysis in the string field theory The pressure components in all p spatial 

directions are respectively: pjx = q (2 < j < p) and 

PlK = PjK + j-^^^PK- (6.4) 

Therefore, pjx — —pK and piK ~ pjK + Pk near the end of the condensation. And so 
PjK — > — oo and pix — > as t ^ oo, which further give wjk — > —1 and wik 0. This is 
consistent with the appearance of a D(p — l)-brane in a vacuum. 

Second, we consider the approach to vacuum of a homogeneous field with an expanded 
inhomogeneous part. The energy density pv and the pressure piv near the vacuum can 
be analysed using the results in Section 3. The final values of pv and piv at the end of 
the condensation are not obviously determined because the two quantities V{T) and 1 + y 
both tend to zero as t ^ oo. They can be estimated in terms of the expression of 1 + y 
in Eq. ( 3.1^ ), from which we see that 1 + y ~ Ae"^^^"'"^^* — > Ae~'^^^ when 6 < (e.g., 



near peaks) and 1 + y — > — CAx^e"^^* when 6 > (e.g., near troughs). To keep I + y 
positive, we assume that C = bQ — bQis negative. For the inverse hyperbolic potential (^), 
V{T) ~ 2Vm exp [—I3{D + t)]. Thus the energy density near the vacuum is 



Pv 



^e-/^^, (6 < 0), 



Thus the energy density tends to constant near peaks and diverges at troughs Ax = as 
t oo. When b is not too large, pv ~ l/|Ax| near but away from the trough at Ax = 0. 

In Fig. ^ we present the plots of the energy density pv and its inverse l/pv relative 
to the plots of the second derivative b = T" at a time t = 16. The figure indicates strong 
relation between py and T": the energy density peaks where T" also peaks. By comparison 
with Fig. I, we see that the sharpest peaks in T" occur at troughs in T. The top plot shows 
evidence that around the peaks of the energy density, we have the approximate relation 
l/pv ~ |Ax|. 

It is interesting to note that the tachyon fluid is not quite pressureless in general. 
Indeed, from Eq. ( |6.2| ), noting that g — > near the vacuum as i — > oo, we find pv'- 
PiV ~ Pv iT)'^ / [I + (VT)^]. Given that T "freezes", we see that the small fluctuations 
give rise to a small pressure. The equation of state parameter Wiv is non-zero. 



7. The generalised DBI action 

The problem with the divergence in the gradient of the field near kinks has been recognised 
for some time ^, IC], and regularisations of the action proposed [^, 35 1 with respect to 
this problem, which have the important property that their static kinks have finite field 
gradients as the minimum potential Vq ^ 0. In what follows, we will consider the modified 
action given in pl in the 1 + 1 dimensions: 



-y(r)(i + y)2 



l(i+.) 



(7.1) 
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Figure 4: The correlation between the energy density pv{t,x) near the vacuum and the second 
derivative T". The plots are made at t = 16. It is clear to see that pv{t,x) peaks at positions 
where T" also peaks and that l/pv ^ |Aa;| near the peaks of pv- 



where e is a small positive value. The equation of motion from it is: 

f = /{2(1 - e^)fV,TV'f - ^(1 + y)(l - ey) 

+(1 + e)[(l + y)V^T - (1 - e)ViTVjTV'V^T]}, 



(7.2) 



where / = !/ (1 + e)(l + ViTV^T - et^) 

First, We consider the field expansion around the kinks and antikinks in the case with 
the hyperbolic potential (^I^). Doing the field expansion like ( |3.lD , we get the equation of 
a{t): 

.. e/3V , (1 - e)(2aa^ + c) , 

« = TT7-TT7 + TT^^ ^^''^ 

So when a grows large and satisfies 1 <^ a <^ l/v^i we have a similar solution to Eq. 
( |3.6| ). When a » the term becomes important and a begins to execute damped 

oscillations, a settles to a = l/\/e. 

To show this point clearly, we can transform Eq. (^^) into a soluble form, assuming 
1 + ~ a^, and neglecting c. When e 7^ 1/2, we set a = z~^/^^~'^^^ and then have: 



z + 



dU{z) 
dz 



where U{z) is: 



U{z) 



e) /32z2 



+ 



1 



0, 
1 

2 ' 



For e > 0, the minimum value of happens at 

1 

a = 



(7.4) 
(7.5) 

(7.6) 
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This value should be the final field gradient T' ~ a as t ^ oo. When e = 1/2, we can set 
a = and get the same equation as in ( [7 .41 ) but with a different potential U{z): 



U{z) = ^ (e^^ - 4z) . (7.7) 

The position according for the minimum potential U{z) is z = or a = \/2, which is 

consistent to the above result of e 7^ 1/2. 

For the expansion around extrema, we get: 



1 + e 



1 - eT^ 



rp2 

° (7.8) 



When e < 1, 1 — eTg is always positive if |To| ^ 1- Similar to the analysis for the DBI 
action with e = 0, |To| still goes to 1 here at the end of the condensation. 

With the approximation Tq ~ 1 for positive Tq, the equation for b{t) can be expressed 

as: 

(1 - e)dt{b - b^) = -2(/? - eb){b - b^). (7.9) 

For e < 1, it can be shown both analytically and numerically that there are still two possible 
final values for b{t) at t — > 00. When b tends to zero, we can write the solution similar to 
( |3.14 ). So it can evolve to zero at t ^ 00. When b{t — > 00) 00, b grows to infinity more 



quickly than the case with e = 0. Thus caustics are possible to form in a finite time for 
the modified DBI action. 

Hence, although this modified effective action solves the problem of diverging gradients 
near kinks, it still has unstable solutions and breaks down as the vacuum is approached. 

8. Conclusions 

We have investigated the inhomogeneous tachyon condensation process, focusing on the 
solutions near kinks (where T = 0) and also in regions where the field is approximately 
uniform. 

We obtained an approximate space-time dependent solution near the kinks and anti- 
kinks, verifying that the spatial derivative of the field T' tends to infinity in a finite time. 
We show that this singularity can be avoided by modifying the effective action in 1 -|- 1 
dimensions, which is consistent with the result obtained in the static case by 35 [. 



We then studied inhomogeneous tachyon condensation in the absence of kinks. In 



(3.18), an analysis assuming that the tachyon field obeyed the eikonal equation T + 
VT^ = 1 showed that near troughs T" diverges in finite time, which was interpreted as the 
production of caustics in free-streaming matter. Directly from the semi-analytical solutions, 
we learnt how the field obeys this equation and how the solutions from this equation lead 
to singularities in T" . A linear analysis shows that perturbations "freeze", except near 
extrema; Adding in higher order terms shows that the curvative (V^T) increases near 
troughs (extrema with V'^T > 0) and decreases near peaks (V^T < 0). For the former 
case, this leads to discontinious gradients near troughs, which in the eikonal approximation 
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leads to caustics. Moreover, the energy density pv diverges as t — >■ oo near troughs. 
The pressure of tachyon matter is smaU, but the equation of state parameter w does not 
generically vanish. 

The same analysis applies also to BSFT effective theories, and to another simple pro- 
posal for modifying the effective action |5|, |3^ . The large gradients and higher derivatives 
signal a breakdown of the effective action ( |1.1D for describing the dynamics of the tachyon 
field near the vacuum. Whether or not caustics actually form remains open for further 
discussion. 
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